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Details of recent developments of the Feynman diagram analyzer DIANA (Diagram ANAlyser) are presented. 
Apart from some discussion about QGRAF and new plotting features, we concentrate on a new sophisticated 
mechanism of momenta mapping. 



The C-program DIANA (Diagram ANAlyser) 
for the automatic Feynman diagram evaluation 
was first documented in pQ 1 . Meanwhile there 
exists a growing series of applications 2 for which 
DIANA was used. 

The current DIANA version is 2.37. Here we 
describe some details of the recent development. 

1. Compatibility with current QGRAF 
versions 

DIANA uses QGRAF |3] as diagram genera- 
tor. By now there exist several QGRAF versions 
which, however, have incompatible syntax of their 
input files. In most cases the user does not need or 
does not want to know which version of QGRAF 
is in use. DIANA takes care of this. It automati- 
cally investigates which version of QGRAF is im- 
plemented. Starting from version 2.37 it also pro- 
vides the user with information about the version 
used for the diagram generation. At present, all 
QGRAF versions 2. • • • and 3. • • • are supported. 
From the DIANA point of view there is no differ- 
ence between them; the differences concern only 
syntax. 

2. New plotting features 

QGRAF cannot generate two-line vertices. 
This causes problems when DIANA is requested 

*On leave from BLTP, Joint Institute for Nuclear Re- 
search, Dubna, Russia 

1 See also 

:/ /www. physik. uni-bielefeld.de/tentukov/diana. html 



to introduce counterterms. In this case it pro- 
duces three-line vertices but does not show one of 
them, i.e. in this case DIANA skips the image of 
the "spurious" line in the diagram. Another new 
feature is introduced to suppress labels. In ver- 
sion 2.35 of DIANA and later ones these options 
(changing the model extensions 0]) are taken into 
account in terms of "optional extensions" in the 
propagators of the model file, e.g. in 

[le,Le;l; FF(num,f num,vec , mle)*i_; 

mmle; nothing, , ; ] 

the last two extensions ; nothing ,0,0;] (sepa- 
rated by ;) have the following meaning: 

1. the new line type of the propagator "noth- 
ing" means that this line will not be drawn. The 
two following parameters in general stand for line 
characterizations and must be included in any 
case. 

2. the second "extension" apparently is empty, 
which means that no label will be drawn - in pre- 
vious versions this meant the "default image" (in 
the above example: le) 

3. Automation of momenta distribution 

Again, starting from version 2.25, momenta can 
be introduced automatically. There are several 
different approaches |5l(i| . 

As "standard" way of attributing momenta, we 
consider the possibility to calculate, e.g., from ex- 
ternal momenta the so called "chords" , character- 
izing the virtual lines apart from the integration 
momenta. 
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A more specific problem is the following: oc- 
casionally topologies have to be generated from 
more complicated ones by scratching lines and 
one may want to stick to the momenta introduced 
for the lines which are kept. This option was de- 
scribed in [5]. 

A new and more difficult situation occurs when 
the user wants to "scratch" some external lines. 
Consider e.g. the Anomalous Magnetic Moment; 
this is a 3-point function of the type fermion - 
fcrmion with an external photon of zero momen- 
tum. Performing a necessary differentiation and 
putting the photon momentum to zero, the re- 
sulting graph appears to be a 2-point function 
with higher powers of scalar propagators. In this 
case we need to map the momenta of the 3-point 
function to the corresponding momenta of a 2- 
point function, putting one external momentum 
to zero. 

At present, DIANA is able to map momenta 
from n-point functions to m-point functions (m < 
n). Here we provide some details. 

Let us consider one of the diagrams relevant for 
g — 2 in two loop approximation, Fig^ 




Figure 1. A typical diagram for the two- loop 
Anomalous Magnetic Moment. Arrows indicate 
the direction of momenta. 



After differentiation w.r.t. to the photon mo- 
mentum <7, we have to put q = and evaluate the 
resulting propagator with the only one momen- 
tum p. 




Figure 2. Momenta decomposition. The initial 
momentum flow is a "sum" of a three-point mo- 
mentum flow a) and a two-point momentum flow 
b). 



To perform the differentiation w.r.t. q we for- 
mally consider at first a diagram with the only 
nonvanishing momentum q, all others zero. In 
a FORM program the differentiation is then per- 
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formed w.r.t. q and after that for zero momentum 
q, momenta fitting to the three-point function are 
introduced, see Fig[21 

In a first step we generate the FORM input 
with momenta as in Fig|5^). After differentia- 
tion w.r.t. q, we look (e.g. in another program 
like MINCER [7j) for the corresponding two-point 
topology (FigEh), add for technical reasons all 
momenta from the two-point topology |2p) to the 
current three-point topology [5^) , put q = and 
generate the FORM input for further evaluation. 

The main problem here is to determine which 
two-point topology corresponds to the current 
three-point topology. In DIANA topologies are 
represented pp in terms of ordered pairs of num- 
bers like (fromvertex, tovertex) (see Fig.[3J). 
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Figure 3. Topology of [2b) as a set of ordered 
pairs of numbers like (fromvertex, tovertex). 

All external legs have negative numbers. The 
momenta distribution according to Fig. [31 e.g., is 
added like 

topology A = 

(-2, 2) (-1,1) (1,3) (4,1) (2 ,4) (3 ,2) (3 ,4) : 
kl+p , kl , k2 , k2+p , kl-k2 ; 

This fixes directions and values of all momenta 
on the internal lines. External lines must be 
known from the process definition. 



A special TM 2 operator was introduced: 

\amputateTopology (varname , tablename) 

It tries to find an amputated topology for the 
diagram under consideration (e.g. [5^, taking off 
q) in the table "tablename" (e.g. from MINCER) 
and returns this topology (e.g. It adds all 

momenta from the table topology to momenta of 
the current topology and saves the resulting 
momenta to the local variable <varname>_M. 

On failure this TM operator returns an empty 
string. If it can't load the table it halts the pro- 
gram. 

Apart from <varname>_M, the operator creates 
several other variables: 

• <varname>_IN - number of internal lines in 
the amputated topology; 

• <varname>_N - name of the found topology. 

The following arrays amputated(original) 

are indices of lines of the amputated topology 
which corresponds to the original: 

• <varname>_E- external lines of original <-> 
external lines of amputated; 

• <varname>_IE internal lines of original <-> 
internal lines of amputated; 

• <varname>_I internal lines of original <-> 
external lines of amputated. 

The absolute value of elements of these arrays are 
indices of the lines and their signs correspond to 
their direction. 

External lines to be removed by the operator 
\amputateTopology are marked by another TM 
operator, 

\rmExtLeg(n) 

where n is the index of the external line. There 
can be more than only one external line which the 
user wants to be removed. 

2 TM language is used by DIANA, see 
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